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Abstract 

Theoretical analysis is presented on quantum state evolution of polarization light waves at frequencies uJo 
and uje in a periodically poled nonlinear crystal (PPNC). It is shown that the variances of all the four Stokes 
parameters can be squeezed. 
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1 Introduction 

During the last decade much attention has been paid to the realization of quantum information protocols [1] 
such as quantum teleportation, quantum cryptography. These protocols are based on the methods of nonlinear 
quantum optics. Nonlinear optical sources [2] play an important role in the generation of nonclassical states 
of light [3] and realization of optical quantum information protocols. The nonlinear optical processes such as 
degenerate parametric down conversion (type I and II processes) [2,3] and the Kerr effect [2,3] are used to 
create nonclassical states of light (squeezed states, polarization squeezed states and entangled states). The 
variance of one of the four Stokes parameters of polarization squeezed light is less than the corresponding 
value for the coherent state. Traditionally the degenerate parametric process (type II) and the Kerr effect 
are responsible for the generation of polarization squeezed states or polarization entangled states in ordinary 
nonlinear crystals. One can achieve suppression of variances of at least one of the Stokes parameters [5-9] 
(5*0, Si, &, Ss) in the type II process by using ordinary nonlinear crystal. Most of the quantum information 
protocols are based on type II process [1] and Kerr effect, which are used to generate entangled states. The 
entangled states are used in the realization of quantum teleporation and quantum cryptography protocols. 
Experiments on quantum teleportation and quantum cryptography are performed by using ordinary nonlinear 
optical crystals with second and third order nonlinear susceptibilities. In the past few years, some experiments 
in the creation and realization of nonclassical and entangled states are performed by using PPNCs with second 
order nonlinear susceptibilities. The PPNCs J10,ll], which have many interesting advantages as compared to 
ordinary nonlinear crystals were proposed by Bloembergen and co-authors in 1962. The main advantages of 
PPNCs against ordinary nonlinear crystals are: the quasi-phase-matching condition between the interacting 
waves; the highest nonlinear susceptibility coefficient can be used; multi-mode interaction of optical waves. 

Recent experiments on quantum noise reduction J12,13] and the generation of entangled states J14] promising 
the applications of PPNCs in the realization of optical quantum information protocols. These experiments were 
based on type I and II processes. It should be noted that the parametric down conversion (type I) and frequency 
sum generation processes have been studied theoretically Jsee for instance, J4] and the references therein] and 
experimentally [12-15] very well. In a PPNC, the possibilities of type II process is much more complicated as 
compared to ordinary nonlinear crystals. The type II process in PPNC can be accompanied by three other 
nonlinear processes J4]. All these nonlinear processes can be quasi-phase-matched at certain coherent lengths 
]4]. Here we will study the generation of polarization squeezed states based on type II processes in PPNC with 
second order nonlinear susceptibility. 

The main goal of this work is to show that PPNCs can suppress all the four variances of Stokes parameters 
below the standard quantum limit. The theoretical work that we present here is the first, to the best of our 
knowledge, to propose PPNCs for the generation of polarization squeezed light. 

The structure of the paper is as follows. Section 2 describes the optical nonlinear processes and their 
Heisenberg equations of motions. Section 3 studies the behaviour of mean photon numbers of degenerate 
polarization (orthogonal) modes at frequencies uJo and tOe. Section 4 analyzes the variances of Stokes parameters. 
The final section summarizes the results obtained in sections 3 and 4. 



2 Equations of motions 



We consider the five-frequency interaction of co-propagating light waves in a PPNC (see Fig.l). The four 
interaction processes of light waves at frequencies u>o, uje, SoJe, SciJe, and 4cJo are [4] 

Ski = fee — kio — fcie + rriiGi = A/ci + miGi, (1) 

UJo + ZoJe = 3uJe , 

5k2 = fcse — fclo — k2e + m2G2 = Afo + m2G2, (2) 

CJe + 3tJe = 4aJo, 

(5A;3 = fee - fcle - fee + JTlgGs = Afe + mgGs, (3) 

2[^e + 2a;e ~ 4[^o, 

5fe = fee — 2fc2e + niiGi = Afe + miGi, (4) 

where Afe=i,2,3,4 is a phase mismatch for an ordinary nonlinear bulk crystal; fen is a wave number of interacting 
modes (n = o stands for ordinary and n — e for extraordinary); Gj = 27rA~^ is the modulus of the vector of a 
reciprocal lattice with period A^; mj = ±1, ±3, . . ., is the quasi-phase-matched order. 

The nonlinear process (1) describes the splitting of a photon of frequency 2a;e into two photons of orthogonal 
polarizations with degenerate frequencies ujo and oje. The second process (2) describes the sum frequency 
generation, i.e. photon of frequency 2ci;e combines with the photon of frequency u)o, which gives rise to a photon 
of frequency 3ti;e. The third (3) and fourth (4) processes are responsible for the generation of the fourth harmonic 
with two different ways, i.e. a photon of frequency lOc and a photon of frequency 3a;e combine, and a photon 
with frequency Aujo appears as a result or the same can be achieved by the combination of two photons with 
frequencies 2aje. 

It has been shown that the nonlinear processes (1-4) can be simultaneously quasi-phase-matched [4] in a 
single domain structure (Gi = G2 = G3 = G4) or at certain coherent lengths I/^oh- The processes (1-4) under 
study can be simultaneously quasi-phase-matched with the following relationship 

Lll = 9iel- (5) 
The processes (1-4) can be described by the following interaction Hamiltonian 

Hi{z) = %(z)[5iaioaieaJj.e'^''"^'' + C2aioa2eaJee'^'°^'' + iso-ieaseCii^e'^''^' + C4a2ea4oe'^''*" + HC], (6) 

where h is Planck's constant, ajn(a^„) is the annihilation (creation) operator of photon of jnth mode at frequency 
jti^n; £,j is the nonlinear coupling coefficient; g(z) is the periodic function equal to +1 or —1 at the domain 
thickness I — A/2; HC denotes Hermitian conjugate. The operators aj„{a'^„) obey the following commutation 
rules 

[aj„,apn] = S.jn,ijn, with j,p= 1,2,3,4 (7) 

The interaction Hamiltonian (6) can be averaged over the period A, if the interaction length z is much more 
than the period of modulation A, i.e. 2 ^ A. Then the interaction Hamiltonian (6) takes the form 

Hi = h[-yiaioaiea2e + 72aioa2ea3e + Isaieaseai + 74a2e"4o + HC], (8) 

where 

+A/2 

7j = ^ J g{z) exp {±iAkjz) = 2^j/{TTmj). 

-A/2 

The Heisenberg operator equations corresponding to the interaction Hamiltonian (8) are given by 

.daio „ + ^ , ^ + _ 

i—, = 72ffl2e"3e "I" 7l "leCS2e , 

dz 

.date , + , , ^ + ^ 

t—. — = 73a3e«4o + 7iaioa2e, 
dz 

.dd2e „ . . , . 

^—, = 27404002^ + 72«lo"3e + Jiaiettlo, 

dz 

= 72aioa2e + 73fflleQ-4o, 

■ ddio ^2 , - - f„s 

= 74a2e + 73aieCl3e- (9) 



Assuming the pump modes at frequencies a;2e, i^4o are classical and non-depleted at the input of a PPNC, i.e. 



a2e ~ A2e, 



(10) 

(11) 



where A2eAo = e"^''^ are the complex amplitudes of the pump modes, we obtain the following linear system of 
equations, after the substitution of quantities (10) and (11) into (9): 



dz 



= -72a3e +7iaie, 



dz 

da^e . , 

-J^ = 72010 + 73«le- 

For simplicity, we introduce a normalized interaction length parameter 

C = ^71- 



(12) 



(13) 



The quantity (13) is introduced into the set of equations (12), which reduce after straightforward algebra to the 
set of equations 



~d^ 



-k2aze + Hie, 



= feflio + kidf^. 



(14) 



where ki = 73/71 and = 72/71. The set of linear equations (14) is solved by applying the Laplace transfor- 
mation: 



aio(C) = All(C)aio -I- Xl2{0^1e + Ai3(C)a3e, 
aie(C) = A2l(C)ai, + A22(C)aie + A23(C)a3"e, 
a3e(C) = A3l(C)aio -I- A32(C)ai"e + A33(C)a3e. 



(15) 



where (0) and 



q = sjl + kl-k 



Aii(C) = -^{—ki coshgC + q cosh + k\) 

q^ 

Ai2(C) = -^(-kik2 coshgC + kik2 + qsinhqO 

g 

Ai3(C) = -i;{-qk2 sinhgi^ + ki coshgC - ki) 

q^ 

A21 (C) = A- (g sinh q( - kik2 + kik2 cosh q() 



A22 (C) = ^ (^2 cosh (jC -f g coshi7C-fc2) 



A23(C) 



;(qki sinhgi^ — k2 coshgi^ + fc2~ 



A31 (C) = '^{1^2 sinh q(^ — ki + ki cosh q() 
A32 (C) = -oilki sinh q( + k2 cosh qC, — k2) 



A33(C) 



(l + g cosh — cosh g(^). 



(16) 



If fcl = fc2 = 0, then the solution (15) corresponds to the conventional degenerate parametric down conversion 
process (type II) [3]. Using (15) we can calculate the statistical properties of interacting modes with frequencies 
uJo, i^e, SuJe at normalized interaction length 



3 Evolution of mean photon numbers 



The evolution of mean photon number of degenerate polarization modes at frequencies uJo and uj^ in a PPNC 
are calculated by using (15) for the following initial conditions at the input of the PPNC: the two orthogonal 
modes are in polarized coherent states |qio >, |aie >, ctio.e = |Q!io,e|e"^i°'= , where |aio,ie:| = \/ < A'"io,e(0) > 
and the wave at frequency 3cJe is in the vacuum state |0 > 

< NloAO > = < "lo! < Oile\ < Q\a'to,e{C)a.lo.e{0\^ > i^le > l^lo > (17) 

The expressions for mean photons (17) are 

< Wi.(C) >= XUQ + A?2(C)kie|' + 2Aii(C)Ai2(C) cos (0i„ + ,^1.) + A?i(C)|aio|', (18) 
< iVi.(C) AL(C) + A^3(C) + AL(C)l«ie|' + 2A2i(C)A22(C) cos (0io + 0i,) + Aii(C)|Qio|'. (19) 

It is well known that the parametric down conversion process (type II) depends upon the initial phases (l>io.e of 
the polarized (orthogonal) coherent states. The values of mean photon numbers (18) and (19) depend upon the 
sum of initial phases, i.e. (f>io + 4>ie. Under the condition (/f>io + = 2ns, the mean photon numbers (18) and 
(19) start increasing rapidly with the growth of interaction length ^, under the condition (j)io + (t>ie ~ tts, they 
start decreasing and then monotnically increasing (see Fig. 2) and s — ±1, ±2, .... 

Fig. 2 demonstrates the dependence of mean photon numbers < Nio > and < A'^ie > on the normalized 
interaction length under the condition (^lo + 0ie ~ i". It is seen that the behaviour of mean photon numbers at 
frequencies uJo (curve 1) and cje (curve 2) are quiet different from the evolution of the mean photon number for 
the case of parametric down conversion (type II) in an ordinary nonlinear crystal. The mean photon number 
(curve 3) for the case of an ordinary nonlinear crystal is calculated by putting ki — k2 — into the expressions 
(18) and (19) under the same initial conditions. 

The difference between the evolutions of mean photon numbers (18) (curve 1) and (19) (curve 2) illustrated 
in Fig. 2 is related with the interaction process (2), which is responsible for the sum harmonic generation, 
i.e, photon of frequency uJo combines with the photon of pump frequency 2uje, which gives rise to photon with 
frequency 3uje. This can be easily seen by evaluating the mean photon numbers < NgeiC) > using (15) under 
the same initial conditions 

< iV3e(C) A32(C) + AL(C)|Qle|' + 2A3l(C)A32(C) COS {ct>lo + 01.) + A^i (C) |q1o . (20) 

The curve (4) of Fig. 2 illustrates the growth of the mean photon number < N-^e > as the normalized interaction 
length increases. 

At earlier stages (normalized interaction length ~ -i- 0.4), the mean photon number of frequency tJe 
decreases (see Fig. 2, (curve 2)) and later starts increasing monotonically. The decreasing of mean photon 
number < A'^ie > is connected with the nonlinear process (3), which is responsible for the sum harmonic 
generation at frequency Aujo and is realized by the combination of photons with frequencies uje and 3uje- The 
fourth harmonic generation process (3) is complicated as compared to (1) and (2). So, the interaction process 
(3) starts acting at later stages of interaction as compared to the interaction process (2). 

The calculations of photon variances < ANi^ ,,{Q > of polarization modes with frequencies lOo and uJe have 
shown that they are super-Poissonian [3]. The same can be seen in the ordinary parametric down conversion 
process (type II) in ordinary nonlinear crystals by putting the nonlinear coupling coefficients fci, equal to 
in the expressions (16) and substituting a3(0) = into the solution (15). 

4 Stokes parameters 

The polarization properties of orthogonal modes at frequencies loo and uje can be analyzed by the Stokes 
parameters [8] 

^o(C) = a^o(C)aio(C) + a+(C)aie(C), 
^i(C) = a^o(C)aio(C)-a+(c)aie(C), 

^2(0 = aUO&ieiO+aUOaioiO, 

MO = il&tiO&ioiO - a^o(C)aie(C)l, (21) 

The Stokes parameters (21) obey the commutation relations of the SU(2) algebra [8] 

[5o(C),Si,2,3(C)] = 0; [5i(0,&(C)] =2iS3(C); 
[S2(C),S3(C)] = 2iSi(C); [53(C),Si(C)] = 2iS2(C). (22) 



The Heisenberg uncertainty relation for the Stokes parameters (21) is given by [8] 

< ASfiO X A5f (C) >> I < Sk > I', fc = 1, 2, 3)(j / J / k). 



'S'0,l(C) = PO±(C) +Pl±(C)«3ea3e + P2± (C) {aieaSe + af^aj^} +P3±(C)"le"le + P4± (C) {aioflSe + ^f^ase} 

+P5±{C){aioaie +a'^oO-le} +P6±(C)"loaio, 

S'2,3(C) = j°'Mao(C){aL ±^3^} + '?i(C){aia3e ± aieajj + q2(C){ai^ ± aL} + g3(C){aioa3e ± a^oOse} 

+q4,{0{aioate ± aio"ie} + g5(C){aL ± alt}], 





(23) 
(24) 



where i 



i stand for 52 (C), S3{Q and 
Po±(C) 



A?2(C)±AL(C)±A23(C) 
P1±(C)=AL(C)±A23(C) 
P2±(C) = Al2(C)Al3(C) ± A22(C)A23(C) 

P3±(C) = a?2(c)±aL(C) 

P4±(C) = An(C)Al3(C) ± A2l(C)A23(C) 
P5±(C) = An(C)Al2(C) ± A2l(C)A22(C) 

P6±(c) = A?i(c)±Aii(c) 

go(C) = Ai3(C)A23(C) 
gi(C) =Al3(C)A22(C)+Al2(C)A23(C) 
92(0 = Al2(C)A22(C) 

g3(C) =Ai3(C)A2i(C) + Aii(C)A23(C) 
g4(C) = Ai2(C)A2i(C) + Aii(C)A22(C) 
g5(C) =Aii(C)A2i(C) 



(25) 



Further for simplicity and clearness, we will write Pj (C)(j=o,i,2,3,4,5,6) as pj and (C)(i=o,i,2,3,4,5) as qj. The 
expressions for the variances of the Stokes parameters (23) and (24) become very lengthy. So, we write down the 
expressions of variances for the Stokes parameter Sj (0 under the same initial conditions applied for expressions 
(17)-(20) 



< ASf (0 >=< S|(C) > - < Sj(C) (j = 0, 1,2,3) 



(26) 



The variances of Stokes parameters (26) are normalized by dividing them by the variance of the Stoke's parameter 
S'o(O), i.e. 



< Ag|(C) > 

<Asm>' 



(27) 



where < AS'o(O) >= |aiop + |cnep. For the case of coherent light, the relative variances (27) take values 1. If 
at least one of the relative variances (27) takes value less than 1, the light is said polarization squeezed. The 
relative variances Vj{C,) read 



-bo± +P2± +P5± + (P2± +2po±P3± +P3± +P5±)|ai<i| 



+P3±\aie\ + 2(p2±P4± + 2po±P5± + P3±P5± + P5±P6± ) |Q!1o | Idle | COs(<j!>io + (pie) 
+ 4p3±P5±|Qlo||aie|^ cos {4>io + 4>le) + 2p5±|aio|^|aie|^ COs2((j!>io + 0ie) 
+ (P4± +P5± + 2po±P6± +P6±)l"lo|^ + 2(P5± +P3±P6±)|Q!lo|^|aie|^ 

+4p5±P6±|aio|^|aie| cos (010 + ((>ie) +pl±\aio\'^ 

-(PO±(C) +P3±(C)l"le|^ + 2p5±(C)|aio||aie| COs((/)lo + (/)le) + P6± (C) I ^lo I ^ ) ^] • (28) 



^2,3(0 



[2q^o + '2ql + qi + 2g| ± 292^101^1'' cos (4</)ie) + (g? 



2q2\aie 



+2(gig3 + 2g2g4 + 2g4g5)|a!io||aie| cos ((pio + </)ie) + 4g2g4|Qio||Qie| cos ((j!>io + cftie) ± 4g2g4|Q!io||Q!ie| cos ((j!>io - 3<?!>ie) 
+4g2g5|aio|^|aie|^ cos (2(/>io + 2(j>ie) ± 2(g| + 2q2q5)\otio\'^\aief cos (2(^io - 2(t>ie) ± 4g4g5|aio|^|cne| cos (3<?!>io - <?!>ie) 
±2q'^\aio\'^ cos (4(j!>io) + 4g4g5|cno|^|Qie| cos {(pio + ^le) + (gs + ql + 4g5)|aio|^ + 2g4|aio|^|aie|^ + 2g5 jaio|* 



± e' 



(29) 



Under the conditions (pio + 4>ie. = tt and k\ < k2 expressions (28) and (29) become more favorable for the 
generation of polarization squeezed light and are evaluated for different initial mean photon numbers in polarized 
coherent states. The Fig. 3-5 demonstrate the evolution of expressions (28) and (29) as functions of the 
normalized interaction length ^. Fig. 3 and Fig. 4 illustrate the relative variances of Stokes parameters So, 2 and 
5*1 and simulate the spontaneous parametric down conversion process, i.e. when the photons with frequencies 
ujo and tOe at the input of the PPNC are each having mean single photons |aio,eP = 1 in their orthogonal 
coherent states. Fig 3. shows that the variances of the Stokes parameters 5o,2 start with the sub-Poissonian [3] 
statistics and later > 0.45 become super-Poissonian. The evolution of the variances of the Stokes parameters 
5*2 and S3 are almost the same. So, that is why we are not demonstrating the variance of S3. Moreover, the 
variances of the Stokes parameters 50,2,3 do not differ from the same variances of the Stokes parameters for an 
ordinary nonlinear crystal. The latter one can be seen by putting the values of nonlinear coupling coefficients 
fci, k2 equal to in the expressions (16) and substituting a3e(0) = into (15). 

The variance of the Stoke's parameter 5*1 is more interesting due to its sub-Poissonian statistics, which is 
shown in Fig. 4. It should be noted that the variance of Stoke's parameter Si for parametric down conversion 
process (type II) in an ordinary nonlinear crystal is super-Poissonian. The latter can be checked by putting the 
values of fci = ^2 = in the expressions (16) and substituting a3e(0) = into (15). So, the PPNC can suppress 
variances of the Stokes parameters 50,1,2 under certain initial values of phases (0io -I- 0ie = tt) of polarized 
modes and nonlinear coupling coefficients (fci < ki). 

The expressions (28) and (29) are calculated under the same initial conditions but with different mean 
photons |aio,ep = 10^ in each polarized mode. The Fig. 5 shows the evolution of variances of all the four 
Stokes parameters 50,1,2,3. It is seen that all the variances of Stokes parameters have sub-Poissonian statistics. 
After analyzing the Fig. 3-5, it is clear that the larger the mean photon numbers in degenerate polarized 
modes, the more the squeezing in the variances of all the Stokes parameter's. Moreover, the variances of Stokes 
parameters 52,3 are ~ 0, i.e. the generation of almost Fock states. 

5 Conclusion 

In this paper we have investigated the problem of the generation of polarization squeezed light in PPNC with 
second order nonlinear susceptibility. In the case of classical and non-depleted modes at frequencies 2uJe and 
4loo, the exact solution of the Heisenberg equations of motions is obtained. The solution is exact quantum 
solution showing new possibilities for the generation of polarized squeezed states of light. We have used this 
solution to calculate the mean photon numbers of interacting modes and variances of Stokes parameters. It is 
found that the evolution of mean photon numbers of interacting modes with frequencies tOo and uie differ with 
the evolution of mean photon numbers of the same interacting modes in ordinary nonlinear crystals (degenerate 
parametric process (type II)). The reasons are found, which explain the pecufiarities of degenerate parametric 
down conversion (type II) in PPNCs. 

Also, it is shown that all the four variances of the Stokes parameters can be squeezed (variances of the Stokes 
parameters are smaller than the value of variance of coherent state) simultaneously. Optimal initial conditions 
are found under which the squeezing of variances of the Stokes parameters can be obtained. 

So, PPNCs can be good candidates for the generation of polarization squeezed states of light, which can be 
used for the realization of quantum information protocols. 
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Fig. 1 Sektch of generation of polarization squeezed light in PPNC with second order nonlinear susceptibility. 
The light waves involved are described by their frequencies luo, uie, lijJe, StJe, and 4ljo. 

3.5 1 \ 1 \ \ 1 \ 1 




Fig. 2 Mean photon numbers < N\o{C) > < ^ie(C) > (2), < ^3e(C) > (4) with frequencies ujo, LOe, and 
3uJe, respectively and < Nio,e{C) > {ki = k2 = 0) (3), as functions of the normalized interaction length The 
curves (1-4) are calculated corresponding to the initial photon numbers < Nio,e{0) >= 1 and < A^3e(0) >= 




Fig. 4 Normalized variance Vi{Q of the Stoke's parameter Si{C,). Curve is calculated corresponding to the 
initial photon number < A'^io,e(0) >= 1 and < A'^3e(0) >= 
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Fig. 5 Normalized variances Vo(C) (0) and Vi((") (1), and V2(C) (2), of the Stokes parameters S'o(C), 'S'i(C) 
and 52(C)- Curves (0-2) are calculated corresponding to the initial photon numbers < A''io,e(0) >= 10^ and 
< NseiO) 



